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A short pedagogical introduction to a differential method used to cal- 
culate multi-loop scalar integrals is presented. As an example it is shown 
how to obtain, using the method, large mass expansion of the two loop 
sunrise master integrals. 
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1. Introduction 

Precision measurements have become one of the central issues in present 
particle physics allowing to test Standard Model and its extensions with un- 
precedented accuracy. To confront them with theoretical predictions it is 
necessary to know two (or more) loop radiative corrections to the mea- 
sured physical quantities. Despite enormous effort of the theoretical physics 
community in this field and existence of many partial results a universal 
method beyond first loop was not developed till now. While the planned 
linear colliders running in gigaZ mode will push the experimental accuracy 
even further. 

One of the promising new directions in the field of multi-loop calculation 
is the differential equations method. A differential equations method, based 
on mass derivatives, has been proposed in [1]. In that approach, amplitudes 
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with a single non vanishing mass m are expressed as a suitable integral of 
the corresponding massless amplitudes, which are taken as known. More 
systematic studies making use not only of masses, but also Lorentz invari- 
ants, as independent variables were initiated in [2] and then continued in 
[3]-[ll]. In [3] small and large expansions were obtained for sunrise type 
master integrals. In [4] the same type of expansions was obtained for the 
two loop two point four denominator master integral. Analytical results for 
pseudothreshold [5] and threshold [6] expansions of the sunrise master inte- 
grals were obtained subsequently. Another application of the method is the 
calculation of the massless off shell double box contributing to 7* — > 3-jet 
process. The differential equations were presented in [7] and subsequently 
master integrals calculated for planar [9] and non-planar [10] topologies. In 
[11] and [12] efforts towards getting 4-loop corrections to g — 2 and 2- loop 
corrections to Bhabha scattering correspondingly were started. The power 
of the method, besides being relatively simple and based mostly on alge- 
braic manipulations, is that its mathematical basis was developed long ago. 
As it will be shown in section 3 the master integrals ever satisfy a system 
of linear differential equations. Theory of such systems of differential equa- 
tions was developed in XIX and at the beginning of XX century, so now it 
is a textbook knowledge (see for example [13]). All that helps a lot, making 
possible mathematical rigor without big effort. 

This paper is organized as follows. In the next section some preliminary 
definitions are given and integration by part identities [14] shortly presented. 
In section 3 it is shown how to get differential equations once the master 
integrals are identified. In section 4 it is shown how to obtain large mass 
expansion of the sunrise master integrals. It is the first application of the 
differential equation method to the calculation of the mass expansions. In 
section 5 a short summary is presented. All algebraic manipulations were 
performed using FORM [15]. 

2. Preliminaries 

For a presentation how the differential equation method works in prac- 
tice in a nontrivial case, but in the same time still not requiring large alge- 
braic manipulations, we have chosen the two point two loop sunrise graph. 
A family of scalar integrals associated with that graph is defined by 



A/ 2 2 2 2 

J (27r)"-2 J (2ttY- 



-02, -a3,/3i,/32) = 



(27r)"-2 {kf + mf )°i (A;| -|- m|)°2 ((p - ki - ^2)^ + m|)°3 



(2.1) 
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where (i = 1,2,3 ) are the masses associated with internal hnes, p is 
the external momentum, ki,k2 are loop momenta and a^, (z = 1,2, 3), 
(3j, {j = 1,2) are integer numbers. The integrals are to be performed 
in n-dimcnsional Euclidean space. It implies we have used dimensional 
regularization and have performed Wick rotation. The scale parameter /j, 
associated with dimensional regularization has been set to 1. Final results 
can be easily rewritten in Minkowski space by changing p'^ —p^. 

Not all of the integrals from the class (2.1) are independent. By means 
of a very simple but powerful method of integration by parts identities [14] 
one can find relations between them. In this particular case one uses the 
relations 



{p-kif^{p-k2r 

{kf + m^)"i (A;| + m|)°2 ((p - ki - k2f + m|)"3 

(2.2) 

where i = 1,2, while denotes one of the momenta p, fci or k2- It is 
crucial that in this way one gets a system of linear equations satisfied by 
the integrals (2.1) with a non-homogeneous terms given by integrals with 
lower number (in this case one) of denominators. This general property 
requires that one should start to solve a given problem from calculation 
of the integrals with smallest possible number of denominators ( or if one 
is lucky enough, one can find them in the literature). The integrals to be 
calculated first for the sunrise problem can be expressed just by one integral 

T{n,m^]= -— — ^ = 7 TTC(n) , (2.3) 

V ' / y (27r)"-2 P + m2 (n-2)(n-4) ^ ' ' ^ ' 

where 



C(n) = (2V^)^*""^r(^3--^j and C(4) = 1 . (2.4) 

With help of (2.2) one finds that only four independent integrals exist within 
the family (2.1) [16]. We choose them as 



Fo{n,ml,ml,ml,p^) = ^(n, m?, m^, mi,p^ -1, -1, -1, 0, 0) , 

Fi{n,ml,ml,ml,p'^) = A(n, mf, ml, m|,p^ -2, -1, -1, 0, 0) , 

F2in,mlrnlmlp') = Ain, mj, ml ml p\ -1,-2, -1,0,0) , 

F^{n,mi,mlmlp^) = A(n,mj,mlmlp'^ ,-1,-1,-2,0,0) . (2.5) 
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Thereafter we will not write explicitly the arguments of the functions 
Fi{n, mf, 771-2, mg, p^) = ^i- The independence of the function Fi means that 
none of them can be expressed by a linear combination of the others and 
polynomials of the function T (2.3) with coefficients in the form of a ratio 
of polynomials in , , , . However an obvious relation occurs 



d 



dmf 



1,2,3. (2.6) 



3. How to get differential equations for master integrals 

Having a limited number of integrals to deal with, which we will call 
master integrals, we can write differential equations they obey. It is as sim- 
ple as to differentiate the given integral and then by means of the integration 
by part identities express the result by the master integrals and the known 
function T. To illustrate as it works we write 



2 d ^ 

I d r d''kid"-k2 1 
2^^dpilJ (27r)2«-4 {kl + ml){k^ + ml){{p - ki - + ml) 
d"-kid'"k2 —p^ + p ■ kl + p ■ k2 

(27r)2'^-4 {kf + ml){kl + ml){{p - ki - k2? + mlf ' 



(3.1) 



The last expression is nothing but a linear combination of three integrals 
from the family (2.1), so one can express each of them as a linear combina- 
tion of the master integrals. This gives after a short algebraic calculation 
[3] 



d_ 
dp 

Similarly one finds [3] 



p^-^Fo = {n-^)FQ + mlFi+mlF2 + mlF^ . (3.2) 



d ^ 

/ D{mlml mi,p2) = ^ m^jF,- + , i = 1, 2, 3 (3.3) 

^ j=o 

where explicit form of functions Tj (which can be expressed by the function 
T (2.3)) and Mjj (polynomials of p^, mf , m^, mg ) can be found in [3]. The 
function D is defined by 
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D{ml,rnl,rnl,p^) = yp^ + {mi + m2 + m^)'^^ + (mi + m2 — ms)^ j 
• (p^ + (mi — m2 + ma)^^ (^p^ + (mi — m2 — ma)^^ )(3.4) 

and is equal to zero at the three pseudothresholds and at the threshold of 
the master integrals. 

Differential equations with m.j as an independent variables can be found 
in a similar way. They are presented in the Appendix of [3] or can be 
obtained from the formulae presented there by appropriate permutations of 
the masses rui, i = 1,2, 3. 

The property that the derivative of a given master integral is a lin- 
ear combination of the master integrals plus terms with smaller number 
of denominators is obviously a general property of all possible multi-loop 
scalar integrals. It is valid due to linearity of the integration by parts iden- 
tities and the differentiation operation itself, and also due to the form of 
the integrands, which are ever ratios of polynomials of masses and Lorentz 
invariants. 



4. An application: Large mass expansion of the master integrals 

Let us assume that one of the square of the masses, say m^, is much 
larger then , and | . The general form of the expansion in that 
region can be written [13] as 



i^ = E(-i)"E^"^7i - (4-1) 

aeA k=0 v"^3J 

where A is a finite set of numbers, whose differences arc not equal to an 
integer number. Other master integrals are related to Fq by (2.6). The 
allowed values of a can be found from the system of equations (with m^ as 
an independent variable) itself. One substitutes the Fi in the system with 
its expansions and by examining the coefficients of the highest powers in 
m^ (they have to be equal to zero) one finds allowed values of a's. In this 
particular case 



^=|-1, ri-3, ^n-4| = |r, si, S2, sgj , (4.2) 

where a shorthand notation was introduced for different values of a's. The 
series with an integer power of a is called the regular series, while the other 



6 



USTR2001 PRINTED ON FEBRUARY 7, 2008 



are called singular series. The singular parts are sources of logarithmic 
terms when expanded around n = 4, the value of n we are interested in. 

Not always all parts of the expansion corresponding to allowed values 
of a's are actually present in the solution. That depends on the initial 
conditions and the regularity of the -Fj at = is crucial for their properties 
[3]. We will see that also in the presented below example. It reflects the fact 
that the differential equations can be satisfied by a wider class of functions, 
not only by the master integrals. It means also that usually one has to 
calculate the master integrals for special values of the parameters by other 
means then the differential equations to fix constants of integration. That 
however is always simpler then the general case. 

Having the allowed values of a one can try to calculate coefficients in the 
expansion (4.1). The crucial point is to find the first coefficient in each of 
the series as the others can be found by solving a system of linear algebraic 
equations (in this case system of 3 linear equations) . Two of the coefficients 
are fixed by non-homogeneous terms in the differential equations and read 



C2 



(n-4)2(n-2)2 



/ 2 2 
(mj^ 777-2 



n-2 
I 2 



(4.3) 



^ C^(77) 
° ((77-4)(77-2))2 



{mi) 2 +{mi) 2 



(4.4) 



The other two cannot be fixed this way as in the non-homogeneous part 
of the equations there is no term ~ (7773)*^ or ~ (7773)*^. One can find 
however the following relations between the next to leading and the leading 
terms in the expansion 



-(77-3) 



{m\ + ml+p^) 

77 



771S2 
^0 



+ -(" - 3)777i(777i -|-p2)-^2 
77 a777f 



d 

d 



= (77 - 3) 
Ff^ = (77-3) 



{p^ -ml + ml) 



d 



dm\ 
d 



rpS2 



^0 



{p' + mi-mi)-^F^^-F^^ 



(4.5) 



Using differential equations (3.2) and (3.3), the expression (4.1) and the 
above relations one finds that Fq ^ satisfies the following system of differential 
equations 
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^ 2 ^ d d 

Prom this system one can deduce, eliminating mf derivatives, that Fq^ 
fulfills a very simple differential equation 



Its solution has the following form 



Fo" = "^'"'"';'"'^' (p-)-°?^+5.(n.n.;.mi) . (4.9) 

where Si(n^m\^m^^ i = 1,2 are still unknown functions. Using (4.7) and 
(4.9) one finds 



d d 

Si{n,ml,ml) = —^Si{n,ml,ml) = , (4.10) 



so the function Si does not depend on masses: S'i(n, m^, 77X2) = 5i(n). This 
information together with (4.5) and (4.6) gives S2{n,m1,m2) = 0. It means 
that Fq^ is a function of n only. Its value will be found later on. 

Similar, but even simpler, analysis can be done for a = S3. From the 
differential equations with as an independent variable one finds that Fq^ 
does not depend on the masses. Using that information, from differential 
equations (3.2) and (3.3) one finds 



2 d 1 



which gives 
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Fo^^=5(n)(/)-^ , (4.12) 

where S{n) is a function depending only on n. As the master integrals for 
n = 4 are analytic functions at = 0, the function S{n) and consequently 
Fq^ have to be identically equal to zero. As all higher order coefficients 
F^^ ~ Fq^, i = 1,- ■ ■, the whole series with a = S3 vanishes. 

The only unknown function Fq ^ (n) can be found using known analytical 
result for Fo(n, 0, 0, m^,p^) [3]. Performing its expansion for large m? and 
comparing the appropriate terms with (4.1), where two masses were set to 
zero, one finds 



F^^n) = 



C\n) 
16 



+ 



(n - 4)2 (n - 4) 



- U+C2 +0(n-4) 



.(4.13) 



As Fo{n,0,0,m?,p^) 



was given [3] in a form of an expansion around 
n = 4 the function Fq ^ (n) is given only in the form of the n = 4 expansion. 
As this was the last missing part of the expansion one can find now the 
complete formula. We report here only leading terms of the expansion, but 
higher order terms can be easily found by algebraic means, if necessary. 



Fo(n, mf, 777.2, m|,p^) — P^^^ terms in (n — 4) 
+C'h{ - ^milog2(mi) + ±mllog{ml) - ^miQ + C(2)) 

+^{rnl + mij log2(mi) 

+ — ^mf + + — 2mi log(mj) — 2ml ^og{ml)^ log(m|) 

( - 3p^ + IGmf log(mi) + 16m| log(m2) - Amf log^(mf ) 
128 \ 

-4ml log^("i2) + (8C(2) - 20)(m? + ml)^ + O ^ - 4^ | (4.14) 

The pole terms being identical to the exact pole terms [3] are not re- 
ported here. As a cross check we have compared the above result with the 
small expansion of Fo(n, mf , m^, ml,^^). We have expanded the first 
two coefficients of the expansion in p^ [3], dependent on mi, m2, ma, assum- 
ing ml ^ ml, ml and found complete agreement between the two results. 
The result (4.14), which is valid for arbitrary p^,mi and ml provided they 
are much smaller then m^, cannot be however deduced from the small p'^ 
expansion itself. 
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5. Summary 

A short introduction to the differential equations method used in calcu- 
lation of the scalar multi-loop integrals was presented. A nontrivial large 
mass expansion of the master two-loop sunrise integrals was obtained almost 
completely by algebraic means. The only 'difficult' task, besides solving 
systems of linear algebraic equations, was to solve two simple differential 
equations (4.8) and (4.11). It shows that the method is extremely powerful 
and opens new possibilities in the field of multi-loop calculations. 
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